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1 Introduction 

In the summer of 1996 LEP started to collect data in the new centre-of-mass en- 
ergy zone corresponding to the VT-pair production threshold and above. During 
the preparation to the first phase of LEP operation |]l[, at centre-of-mass energies 
comparable to the Z mass, it was advocated that different classes of radiative cor- 
rections may turn out to be essential in the interpretation of experimental results. 
Even though it is generally accepted that, if possible, a purely analytical approach 
is more convenient in comparisons of theoretical predictions with data, in many 
cases such as a r polarization measurement 0, Monte Carlo (MC) modelling of 
the observables, including the complicated interrelation of experimental cuts and 
strongly peaked multiphoton phasespace can be handled only with the help of a 
high-precision MC simulation. In the case of the luminosity measurement at LEPl 
[0], QED corrections constitute, even today, the systematic uncertainty surpassing 
the experimental error and limiting the physical significance of the measurement of, 
for instance, the number of neutrino species. 

It is thus of high practical importance to ensure that a similar situation will not 
repeat itself at LEP2 and, if necessary, to take appropriate steps in advance. For 
this purpose a whole family of MC programs and semi-analytical calculations are 
being developed. A comprehensive report on the status of this work can be found 
in the report of the LEP2 Workshop @]. 

Ultimately in PF-pair production and decay we expect to reach an experimental 
precision of order 1%.Q To ensure that uncertainties of theoretical predictions will 
be negligible, it is thus necessary that they do not surpass 0.5% and preferably are 
limited by a 0.3% threshold. 

The purpose of the present paper is therefore to propose a working scheme of 
estimating these uncertainties. Specifically, we want to supplement technical preci- 
sion tests of the MC programs with a discussion of the physical precision. For the 
KORALW version 1.02 it is performed at the technical precision level of better 
than 0.1% (in fact better than 0.01%) for the VT-resonant diagrams. Similar tests 
for the MC program grc4f [^ were performed at the 0.2% level 0] on the total cross 
sections for the W^-resonant diagrams. This technical precision can be improved just 
by spending a lot of computing time for the simulation as it comes from the statistics 
of generating events for the comparison. This way, having specified both technical 
and physical precision of the MC calculation, one will be able to establish the total 
theoretical error of the results. 

We will discuss the general strategy of estimating the physical precision as well as 
illustrate it on the following observables, all in the sscc decay channel: the total four- 
quark (four-jet) cross section, two-quark/jet invariant-mass distribution and cross 
section in the case when the other two may escape detection. We will discuss in 
turn effects such as Coulomb correction, running W and Z widths, effective coupling 



"'^One has to remember here that 'universal' experimental precision does not exist. It depends 
on the measured quantities. 



constants, etc. We will also comment on the uncertainty related to the remaining 
classes of electroweak (including QED) corrections, but we will leave aside questions 
related to QCD final-state corrections and hadronization phenomena. 

Since the publication of Refs. j^, |^ which were based on the version 1.03 of 
the KORALW code, the KORALW program has undergone a series of upgrades, resulting 
in the recent version 1.210. The largest improvements are two sets of complete 
tree level amplitudes for all the possible massive CC-type final states generated by 
GRACE packages versions 1 and 2 [§], interfaced to the KORALW generator as optional 
external matrix elements (amplitudes of GRACE v. 2 are similar to those implemented 
in the grc4f ), along with a new internal presampler to facilitate the general four- 
fermion phase-space generation. For cross-checks the CC-03 differential distributions 
and total cross sections of Ref. were reproduced with the help of this new 
presampler of KORALW. This technical comparison with semi-analytical integration 
gave an agreement at the level of 0.03% at least. 

The grc4f has a kinematical mapping routine, which can treat the singular be- 
haviour of four-fermion amplitudes, and the numerical integration is done by BASES 
[^]; on the other hand, the unweighted event generation is performed by SPRING 
[^]. Therefore by comparing the results from these two packages, the technical test 
of the independent part of the two packages, i.e. the numerical integration part, 
was performed. The gauge-parameter independence of the GRACE amplitudes has 
been confirmed at randomly selected phase points up to the quadruple precision of 
the variables. Further comparison was carried out with the results of the algebraic 



manipulation package CompHEP |Ty, |Tl| for the various four-fermion final states and 
a precision of at least 8 digits was found for series of randomly chosen points in the 
phase space. 

To cross-check the KORALW and grc4f programs, we compare in Table 1 the 
numerical results for total cross-section for various input parameter configurations, 
as described in the table caption and later in the text. This constitutes a significant 
test of the two programs, as agreement is whenever expected better than 0.3%. The 
remaining numerical results of the paper are produced by KORALW and cross-checked 
by grc4f wherever possible. 

Finally, there is one more highlight in this paper. In the course of various four- 
fermion simulations we found an unexpected and strong artificial effect due to the 
interplay of spin correlations and, seemingly simple, cut-offs on four-fermion phase 
space. We will illustrate these effects with the example of fake peaks in two-quark 
invariant-mass distributions. 

The layout of the paper is as follows. In the present Sect. 1 we have already 
discussed the technical tests of the codes. In Sect. 2 we define our observables and 
elaborate on the question of physical precision. In Sect. 3 we discuss and illustrate 
the spin correlations issue. We conclude our paper in Sect. 4. 



^Available from tittp://hpjiiiiady.ifj .cdu.pl/programs/prograins.html 



2 Numerical results and general scheme 
of calculating uncertainty 

To begin with, we want to stress that in principle nothing hke an overall uncertainty 
of the MC program can be defined. For any new observable, one has to repeat anew 
the numerical analysis as presented here. Nonetheless, from the point of view of 
presenting a general scheme of how to calculate the uncertainty of the program, the 
choice of observable is secondary. 

For the purpose of the present discussion we choose the cess final state (CC-43 
type process) of VT-pair decay because of its relatively high cross section. We take 
the following observables: (i) total "visible" ((T4) cross section for the four-fermion 
final state and (ii) the invariant-mass distribution of ss, in the case where cc are 
escaping detection and the corresponding integrated cross section ((T2). We will call a 
fermion "visible" if its transverse momentum is above 10 GeV and | cos 6beam\ < 0.96. 
Otherwise we call it "escaping detection". Our motivation for such a choice is the 
following: (i) it can be realized in practice by most detectors, (ii) it excludes jet-like 
activity in the initial state, such as off-mass-shell initial-state photon bremsstrahlung 
(or initial-state jet activity in the framework of the phenomenology oi pp colliders). 

In general, we will follow the recommendations of LEP2 Workshop 0] for input 
parameters setting, but as it is essential for this paper, we will explain in detail all 
input parameters and switches of KORALW and grc4f in Appendices A and B. 

As a centre-of-mass energy, we take 161, 195, and 350 GeV, corresponding respec- 
tively to VT-pair production threshold, maximum expected centre-of-mass energy of 
LEP2, and ti threshold. 

2.1 General scheme of calculating uncertainty 
and numerical examples 

In the following we will present the general scheme for calculating the theoretical 
uncertainties of KORALW and grc4f MC predictions for observables in final states 
of the T^-pair decay type. We will follow the general scheme that we have already 
applied in other cases 0, |l^, |l^, |14[. For this purpose we need to collect numerical 



results obtained at different levels of approximation. This will enable us to calculate 
the size of different classes of corrections. Approximation levels are listed in the 
caption of the Tables 1,2 and the corresponding detailed information on the KORALW 
and grc4f input is given in Appendices A and B. 

Let us start with the predictions at the lowest approximation level marked as 
input number^ -2-, see Appendices A and B. In this case we use the constant width 
approximation, on mass-shell relation sin^ 9^ = I — M^r/M"^ and tree level doubly 
resonant M^-pair production and decay diagrams (CC-03) only. The results for our 

^ Input number -1- corresponds to the case when the W spin effects are switched off and is 
discussed later on. 



observables are collected for different centre-of-mass energies in Table 2a-c, entry 
-2-, and plotted in thick lines in Figs. 1-6. 

Table 1: Total cross sections [pb] without any cuts for cess production from KORALW 
and grc4f. Different input parameter settings for entries Nos. -2- to -5- are explained in 
Appendices A and B. The energies shown are: 161 GeV (la), 195 GeV (lb) and 350 GeV 
(Ic). Please note numerical differences (~ 0.3%) due to different Coulomb correction 
implementation in KORALW and grc4f, and also due to e.g. YFS formfactor ~ 0.7% 
which is switched on in KORALW. 





la, 


atot [pb], CMS 


energy 161 GeV. 


No. 


KORALW 


grc4f 


Comments 


2. 
3. 
4. 
5. 
6. 


.5750 ± .0002 
.5676 ± .0002 
.6271 ± .0002 
.6194 ±.0002 
.4622 ± .0002 


.5743 ± .0007 
.5660 ± .0009 
.6273 ± .0028 
.6174 ±.0006 
.4565 ± .0005 


CC-03 
CC-03±leading loop corr. 

CC-43 

CC-43±leading loop corr. 

CC-43±leading loop corr. ± ISR 



lb, atot [pb], CMS energy 195 GeV. 



No. 


KORALW 


grc4f 


Comments 


2. 


2.2503 ±.0007 


2.249 ± .002 


CC-03 


3. 


2.0981 ±.0007 


2.095 ± .002 


CC-03±leading loop corr. 


4. 


2.3274 ±.0008 


2.323 ± .003 


CC-43 


5. 


2. 1707 ±.0007 


2.163 ±.002 


CC-43±leading loop corr. 


6. 


1.9753 ±.0009 


1.947 ±.002 


CC-43±leading loop corr. ± ISR 





Ic, atot [pb], CMS 


energy 350 GeV. 


No. 


KORALW 


grc4f 


Comments 


2. 
3. 
4. 
5. 
6. 


1.3771 ±.0005 
1.2771 ±.0005 
1.4132 ±.0006 
1.3101 ±.0005 
1.3616 ±.0006 


1.377 ±.002 
1.273 ±.002 
1.411 ±.002 
1.304 ±.001 
1.344 ±.001 


CC-03 
CC-03±leading loop corr. 

CC-43 

CC-43±leading loop corr. 

CC-43±leading loop corr. ± ISR 



Let us now proceed with different classes of corrections. First let us list those 
classes of higher order corrections that are expected to be enhanced with respect to 
^ ~ 0.2% size. They include: 

1. Coulomb interaction of the H^-pair close to the production threshold, 

2. running Z and W widths. 



Table 2: Total cross section atot [pb] without cuts and cross sections a2 and a^ [pb] of 
two or four "visible" fermions defined as in Section 2. Different input parameter settings 
as explained in Appendix A and B are used for entries No 1 to 6. The energies shown 
are: 161 GeV (2a), 195 GeV (2b) and 350 GeV (2c). 







2a, CMS energy 161 GeV. 






(^tot 


0-2 


0-4 


Comments 


1 


0.57519 ±.0002 


.00041 ± .00000 


0.50269 ± .0002 


CC-03 no spin 


2 


0.57504 ± .0002 


.00035 ± .00001 


0.50785 ± .0002 


CC-03 


3 


0.56762 ± .0002 


.00034 ± .00001 


0.50152 ±.0002 


CC-03±lead. loop corr. 


4 


0.62714 ±.0002 


.00885 ± .00002 


0.52071 ± .0002 


CC-43 


5 


0.61937 ±.0002 


.00877 ±.00002 


0.51429 ± .0002 


CC-43±lead. loop corr. 


6 


0.46219 ± .0002 


.01022 ± .00005 


0.36910 ±.0002 


Like 5 but with ISR 







2b, CMS energy 195 GeV. 






0"tot 


0"2 


0-4 


Comments 


1 


2.25006 ± .0006 


.00480 ± .00003 


1.83470 ±.0006 


CC-03 no spin 


2 


2.25028 ± .0007 


.00382 ± .00003 


1.85607 ±.0007 


CC-03 


3 


2.09812 ±.0007 


.00359 ± .00003 


1.72991 ± .0006 


CC-03±lead. loop corr. 


4 


2.32744 ± .0008 


.00989 ± .00004 


1.90338 ±.0008 


CC-43 


5 


2.17071 ± .0007 


.00950 ± .00004 


1.77361 ± .0007 


CC-43±lead. loop corr. 


6 


1.97533 ±.0009 


.01022 ± .00008 


1.61557 ±.0008 


Like 5 but with ISR 







2c, CMS energy 350 GeV. 






0"tot 


0-2 


0-4 


Comments 


1 


1.37670 ±.0004 


.03055 ± .00007 


0.76319 ±.0003 


CC-03 no spin 


2 


1.37708 ± .0005 


.00713 ± .00004 


0.77531 ± .0004 


CC-03 


3 


1.27712 ±.0005 


.00662 ± .00004 


0.71982 ± .0004 


CC-03±lead. loop corr. 


4 


1.41318 ±.0005 


.00913 ±.00005 


0.79474 ± .0004 


CC-43 


5 


1.31010 ±.0005 


.00853 ± .00005 


0.73714 ± .0004 


CC-43+lead. loop corr. 


6 


1.36161 ±.0006 


.00871 ± .00005 


0.79493 ± .0005 


Like 5 but with ISR 



3. effective coupling constants which affect the relation: sin^^^ = 1 — M^y/M"^, 

4. YFS formfactor and QED Initial-State Radiation (ISR) for the KORALW or 
QEDPS [H for the grc4f . 

It is rather straightforward to include first three types of corrections in the case 
of CC-03 amplitudes and the appropriate results are collected in Table 2a-c under 
entry -3-. In this case we follow a recommendation of the LEP2 Workshop [Q for 
so-called CC-03 diagrams in the choice of definition of sin^ 9w (see Appendix A for 
details) and in switching to s-dependent Z and W widths. As we can see, corrections 
are typically rather small. We thus expect higher-order corrections to be negligible. 



Let us also point out that the dominant higher order corrections of the vacuum 
polarization type are already included in definitions of running widths and effective 
coupling constants. The systematic uncertainty related to those corrections is thus 
included in the uncertainty on the numerical values of the program input parameters 
and not separable from them. 



The Coulomb correction implemented in KORALW is taken from Refs. |jT6|, |T] 
whereas grc4f takes it from Ref. [|T8|. This leads to numerical results that differ by 
about 0.3-0.4% for all energies. According to [Q (vol. 1, p. 117-119) the deviating 



formulae of Refs. ||T^ and |TB[ constitute equally well justified representations of 



the Coulomb phenomenon; however, see also Ref. [0 for a detailed discussion. 
Nonetheless this 0.3-0.4% difference will be reflected (indirectly) in our flnal estimate 
of physical precision. As for the higher-order corrections beyond 0{a) to Coulomb 
correction, it is shown (see e.g. |17, 0) that they are numerically small, below 0.1% 



and can thus be neglected in our considerations. 

Let us now turn to the effects of complete tree level amplitudes. To this end 
we have to switch off all corrections introduced in entry -3- and return to the input 
deflnition as in entry -2-; we then implement the complete set of tree level spin 
amplitudes, the so-called CC-43 graphs. We collect results under entry -4- of Tables 
1 and 2 and in the thin line of Figs. 2, 4 and 6. As we can see, the cross-section 
increases by a few per cent. This is mainly due to the inclusion of the Z-7 interme- 
diate state. Note that when complete Born level spin amplitudes are included, i.e. 
in the presence of a Z peak, another, Z-like peak, in ss invariant-mass distribution 
at higher centre-of-mass energy, looks particularly suggestive, see Figs. 2, 4 and 6. 
The cross section for this faked resonance is 0.0015 pb at 195 GeV, which translates 
to about 2 or 3 events per overall time of operation of every LEP experiment if 
the summation over all hadronic flnal states is performed. (We will return to this 
phenomenon in Sect. ^. 

The simultaneous inclusion of the complete set of tree-level diagrams and incom- 
plete 0{a) corrections, for example running PF-width, as in entry -3-, breaks the 



gauge invariance. See e.g. [19| and also [20, 21, 22|. The other corrections of entry 



-3- may also cause numerical effects upon combining 'by hand' with background 
graphs. Until a complete 0{a) calculation for the four-fermion processes is done, 
the question of coherent extension of partial 0{a) corrections to the background 
case will always be present. 

For the time being we proceed as follows. Formally speaking, in order to combine 
corrections of entries -3- and -4- of our Tables 1 and 2 it is necessary to perform 
several runs of the MC and later to add the corresponding corrections according to 
the following formula: 

X = X2 + (X3 - X2) + (X4 - X2) + h.o.t. (1) 

where Xi correspond to the different approximations (of observable X) as summa- 
rized in the caption of Tables 1 and 2. Equation (|l]) is nothing more than the 
multiparameter linear interpolation based on Taylor expansion in some 'parameter 

6 



space' (symbolically denoted by /'s), in which each direction corresponds to a given 
set of 'orthogonal' corrections (-3- and -4- for example): 

f)X f)X 

AX(/3, U) = —A/3 + —A/4 + . . . = AX3 + AX4 + h.o.t. (2) 

C/3 CIJ4 

and the missing terms are of higher order in the expansion. Approximation (|T]) 
therefore gives a safe way of combining corrections. 

On the other hand this solution may be rather painstaking, as in the case of an 
observable including complicated experimental cuts, several simulation runs for the 
same processes but different theoretical input may be necessary. This may take an 
enormous amount of CPU time if the full detector simulation is requested. 

For this practical reason we propose a different way of combining corrections. 
It simply amounts to switching all of them simultaneously. This common sense 
interpolation is of course Monte Carlo dependent, as each code most likely uses 
different ad hoc method of combining corrections. The ultimate test of common 
sense interpolation can be done by comparison with the linear interpolation of eq. 
(|ID. If these two interpolations agree within the required precision, then one can 
safely use common sense interpolation. However, this has to be checked for each 
type of observable separately, just as we do for sscc ones in here. Specifically, in 
entry -5- of Table 2 we give the result of simultaneous switching on corrections -3- 
and -4-. We can see by direct inspection of Table 2 that our eq. (|l|) reproduces at 
0.3% precision level entry -5- from the corresponding entries -2- and -4-. Strictly 
speaking, this is true for total and a^ cross sections. For o"2 the agreement is at the 
1% level. However, as the a2 is typically two orders of magnitude smaller than the 
others, the 1% agreement is more than sufficient. We can therefore conclude that, 
in this case, our common sense interpolation approach of simultaneous inclusion of 
all corrections works at the required precision level. 

So far, our results were obtained with ISR and YFS form factors switched off. If 
we add these effects, we obtain finally predictions including most of the necessary 
corrections.^ We collect results in entries -6- of Table 2. 

This way, we have exhausted our list of 0{a) enhanced corrections, and we are 
left with the remaining, genuine 0{-) ones. To state it simply, before a complete 
0{a) calculation of corrections to e"*"e~ —>■ 4/ are available, the general uncertainty 
of the C(-), due to QED/EW corrections, must be assumed in the results. 

At this moment in order to quantify the above vague statement we will use a 
number of partial results on 0{-) corrections. We will briefiy cover them in the 
following. 



^ A brief discussion of QED initial-state corrections was presented already at the LEP2 Work- 
shop. One has to keep in mind that even though kinematical configurations with an arbitrary 
number of photons are generated by our programs and corresponding uncertainties were discussed 
for KORALW in Ref. Q, and for grc4f in Ref. |Q, matrix element is limited to the second-order 
leading- logarithmic accuracy in both cases (on top of correct all-order soft limit!). 



1. We briefly reviw the formulation of QEDPS that is used in the grc4f . The 
algorithm is completely in parallel with that of the parton-shower model in 
perturbative QCD, which has been well known for a long time. The basic 
assumption is that the structure function of an electron, with the virtuality 
Q^ and the momentum fraction x, obeys the Altarelli-Parisi equation 

in the leading- log (LL) approximation P3 . This equation can be converted to 



D{x, Q') = ^(Q^ Ql)D{x, Ql) + ^j ^n(g2, K') J ^P{y)D{x/y, K'), 



the integral equation 

"^ TT/^2 r^2N f dy 
y 

(4) 
where 11 is the Sudakov factor. Here, rigorously speaking, Q'j. should be ml 
as it gives the initial condition. For simplicity the fine structure constant a is 
assumed not running with Q^. 

The integral form eq. (0) can be solved by iterating the right-hand side in a 
successive way. Then it is apparent that the emission of n photons corresponds 
to the n-th iteration. The information on the transverse momentum can be 
obtained by solving the kinematical equation. 

In the formulation there appear two parameters, Q^ and Qq. In the program 
the following values are chosen: 

Ql = mle = ml X2.71828---, Q^ = lO'^^ GeVl (5) 

The former value was settled to effectively take into account the constant term 
—1 of (3 in such a way that (3 = (2a/7r)(ln(s/mg) — 1) = (2a/7r) ln{s / (mle)) . 
Since the second parameter is unphysical, any physical observable should not 
depend on it. It has been checked that increasing Ql up to O(mg/10) leaves 
the result unchanged in the statistical error of the event generation. 

We conclude that, in total, the related uncertainty is (a) 0.1% from comparison 
between the QEDPS and the structure function method and (b) 0.65% due to 
neglecting the overall K factor from the exact 0{a) calculation [p^. 

Let us now go back to the YFS scheme and to the discussion of the main un- 
certainty (in that language) related to initial-state bremsstrahlung. The YFS 
formfactor is a consequence of the rigorous and general program of resumma- 



tion of soft photons in any QED process, as given in the classical paper ||25 
In that sense it is universal and process-independent. On the other hand, it 
originates from a particular (although well motivated) choice of approxima- 
tion in extracting 'soft' parts from the Feynman graphs. It can thus be argued 

8 



that there is some arbitrariness in its definition. In s-channel reactions such as 
muon or tau pair production, this formfactor, which numerically constitutes 
+0.7% overall correction, had to be included. This was based on the direct 
calculation of QED corrections up to the second order in a. In the case of 
PF-pair production such calculations are missing so far, and there are opinions 



26| that such correction (YFS formfactor) should perhaps not be included be- 
cause in this case cancellation of terms ~ - ■ tt^ from different sources cannot 
be understood. 



3. In Ref. |27|, the YFS resummation of real and virtual soft photons is extended 
to the case of bremsstrahlung from heavy bosons W. This is done for massive 
W^s with finite width in a manifestly gauge-invariant way. The results of the 
calculation are then implemented in the four-fermion Monte Carlo program 



YFSWW2. The numerical results presented in Ref. |^ show that, for the 
total cross section, the net effect of real+virtual emission from IV-pair, after 
subtracting the Coulomb correction, is of the size of 0.4% for LEP2 energies 
and 0.8% at 500 GeV. The size of this correction, absent in KORALW, should be 
included in the final physical uncertainty. 

4. Another technique of computing ISR in a gauge-invariant way has been de- 
veloped in |^8[. The t-channel neutrino was split into two oppositely fiowing 
charges. These charges were then ascribed to initial and intermediate (ly's) 
states respectively. In this way, certain terms effectively have been rearranged 
between initial- and intermediate-state-type corrections. Numerical results 
presented in ||2^ are: 0.4% for LEP2 energies and 1.5% for 500 GeV. This is 



in qualitative agreement with the results of Ref. [p7| . 
5. The interference effects to the process of VT-pair production and decay have 



been analysed in |£9|, [SQ]- It has been shown therein that, for sufficiently 
inclusive quantities (e.g. total cross section) these effects, both real and virtual, 
are suppressed by an additional Tw/Mw factor with respect to the genuine - 
corrections. 

6. Finally, another possible way of estimating QED/EW uncertainties is to use 
the results of the calculation of the complete 0{a) corrections to stable on- 
shell ly-pair production of Refs. [^, |3^. These corrections are, however, 
to our knowledge, not implemented in a full-scale four-fermion Monte Carlo 
program, which makes direct comparisons more difficult. Based on these on- 
shell results the uncertainty of the four-fermion total cross-section at LEP 2 
has been estimated in Ref. g (vol. 1, p. 127) to be 2% for 161 GeV, 1-2% 
for 175 and 190 GeV. 

In the following we will mostly rely in our estimates of the QED/EW uncertainty 
on the results of Ref. p^. Specifically we take as the estimate the size of the 



YFS form factor for the ISR, summed linearly, as discussed above with the size 



of intermediate state (Pi^'s) YFS real+virtual corrections of [^. On top of that 



we impose a safety factor of 2. This procedure yields 2 x (0.7% + 0.4%) ~ 2% 
and we quote a total QED/EW uncertainty of 2%. (The study of the QEDPS also 
gives 0.7% uncertainty from comparison between the QEDPS and the exact 0{a) 
calculation.) Let us note also, that we did not explicitly include the 0.4% uncertainty 
due to Coulomb correction discrepancy between Refs. |T^ and [jl8[ mentioned earlier. 



We include it in the safety factor of 2 (it is in a sense related to the PF-state 
uncertainty) .0. This is not an improvement with respect to the similar precision 
estimate for the total cross section based on on-shell experience (2% for 161 GeV 
and 1-2% for 175 and 190 GeV) given in Ref. [| (vol. 1, p. 127). 

We hope that a careful study of terms proportional to - -vr^ with special emphasis 
on the matching terms from different calculations of partial 0{a) contributions may 
lead to a reduction of this systematic error. However, the final answer may have to 
wait for results of complete 0{a) calculation. 

Another uncertainty, related to higher-order parts of our ansatz on initial state 
bremsstrahlung matrix element, we take as 0.1% and discard. On the basis of 
inspection of our Tables, we believe that our ansatz on the simultaneous inclusion 
of all corrections even though it violates in principle gauge invariance, in practice 
introduces for our observables an uncertainty of order of 0.1% which we discard as 
well. 

The other component of the total uncertainty to be mentioned here is the ques- 
tion (briefly covered in the introduction) of the technical precision of KORALW and 
grc4f codes. For KORALW it was discussed already in and for grc4f in [^]. Numer- 
ical results of the present paper were reproduced by the two programs and agreement 
was found at the level of 0.3%. Based on all these factors, we conclude that in the 
case of our observables the technical precision of both the KORALW and grc4f pro- 
grams is not worse than 0.3%, and therefore negligible in comparison to the physical 
one. 

We believe that we can sum all the above uncertainties in quadrature, obtaining a 
final result for the uncertainty of K0RALW/grc4f predictions, including technical and 
electroweak effects, for our observables of 2%. The dominant source of uncertainty 
originates from the lack of complete 0{a) calculation of electroweak corrections to 
combined I^-pair production and decay process. 

We have excluded from our considerations an uncertainty in the program nu- 
merical input parameters, such as the Z and W masses and widths as well as the 
W branching ratios. Note that in Pi^-width and W branching ratios there are hid- 
den total final-state QCD corrections. The related uncertainty is not included here 
either. 

We have excluded from our discussion yet another source of systematic error 
which is related to hadronization, jet definition, QCD perturbative effects in the 
final state, etc. Note that this includes final-state QED bremsstrahlung from quarks 

^However a direct summation of all three uncertainties gives 1.5% bottom-rock precision level. 
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which cannot be separated from jet formation. These effects may constitute size- 
able corrections, especially at higher energies. For the purpose of simulation of such 
effects we use the JETSET algorithm and we refer to [^] for a discussion of therelated 
uncertainty. As suggested in Q] (vol. 2, p. 172; "good bet" approach) we leave this 
topic to independent research. We stop at generating a random choice for colour 
recombination between charged-current (CC) and neutral-current (NC) configura- 
tions according to the size of corresponding separate matrix elements squared. So 
far, we neglect interference effects between the two configurations, as well as related 
Bose-Einstein effects. For more details, see Ref. |^ vol. 1 p. 124 and 190. We 
believe, however, that effects and uncertainties related to QCD/hadronization may 
require additional separate studies to be completed in the near future ||33[| . 

3 Spin correlations 

In this section we will show the importance of spin correlations in the presence 
of cut-offs in Pi^-pair production on the example of e^e~ -^ sscc. This is to be 
contrasted with the 'no correlations' case of e~^e~ -^ W^W~ -^ sscc where VT's are 
produced 'on shell' and then independently decayed. To this end we return to the 
'lowest approximation level', as described in Sect. p!.l| , marked as entry -2- in Table 
2 and shown in thick lines in Figs. 1, 3 and 5. We now switch the spin correlations 
off. The results for the same observables and cut-offs as before are shown in Figs. 
1, 3 and 5 (thin lines) and in Table 2 (input number -1- ). 

Let us point again to the spectacular peak in the ss invariant-mass distribution, 
which becomes more and more profound for centre-of-mass energies higher above 
the WW threshold. It is present in the case -2- (thick line) but not (!) in the case 
-1- (thin line) when spin correlations are switched off. It is therefore a genuine effect 
of an interplay of our veto cut-off on c-quarks with spin correlations in W decays! 
This exercise proves that any kind of 'on-shell' approximation of this kind may lead 
not only to quantitative few or several per cent inaccuracies, but, upon applying 
cut-offs, to misleading qualitative changes in the overall picture. 

Note that the cross cection for our faked 'object' is of the order of 0.0015 pb for 
the cess final state alone. This translates to 2 or 3 such events per LEP collaboration 
if all hadronic final states are taken into account. 

Moreover, it should be kept in mind that our choice of observables and cut-offs 
is in a sense random and that similar effects can be expected in other distributions 
as well. 



4 Summary 

In this paper we have discussed different types of non-hadronization corrections for 
the three observables: 'visible four-quark' total cross section and invariant mass 
distribution and cross section of 'visible two quarks' from the sscc final state at 
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LEP2 and higher energies. We have found that in most cases MC simulation can 
incorporate all necessary corrections, such as Coulomb corrections, effective coupling 
constants, running Z and W width and complete set of all tree-level diagrams. Even 
if in principle it can be shown that such a solution leads to gauge non conservation, 
we have found that for practical purposes, as in the cases presented above, all 
corrections can be included simultaneously. This common sense interpolation works 
because of our reasonable/lucky choice of the gauge and observables discussed in 
this paper. This is of practical importance as it eliminates the need to calculate 
corrections separately in several runs of detector simulation. 

We conclude that the overall uncertainty of theoretical predictions for our ob- 
servables is 2% with dominating contribution from the lack of complete 0[ci) calcu- 
lations of combined W^-pair production and decay process. This is already a factor 
of 4 short of the ultimate goal for LEP2 phenomenological preparations, even though 
we still exclude further QCD/hadronization uncertainties. It is, however, sufficiently 
precise for the statistics of the first year of LEP2 operation. 

On the other hand, we have found that seemingly innocent and natural choices 
of cuts (which can be motivated by detector or background elimination needs) may 
lead to a very strong deformation of the 4-fermion signal and even to unexpected 
faked peaks strongly enhanced by spin effects. We expect that a careful study of 
MC predictions, including complete spin effects, may be necessary in many cases 
such as four-jet/leptons final states or six-jet /lepton final states (e.g. involved in ti 
production and decay). This is important not only in the case when a given channel 
is studied by itself, but also if it is treated as the background. 
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Tateya collaboration, J. Fujimoto, T. Kaneko, K. Kato, S. Kawabata, T. Munehisa, 
N. Nakazawa, D. Perret-Gallix, Y. Shimizu and H. Tanaka. They are supported in 
part by the Grant-in-Aid(No. 07044097) of Monbu-sho, Japan. M.S. and Z.W. are 
pleased to thank S. Jadach, W. Placzek and B.F.L. Ward for numerous discussions 
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Appendix A 

In this appendix we carefully review input parameters as used in the production 
of our Tables with the help of KORALW. 

1. * set 1 * 

We will start with parameters as set for entry -1- of our Tables and later on 
we will concentrate on respective differences only. The exact meaning of the 
parameters is described in the KORALW documentation, see Ref. |^. 

500000 NEVTDT 

KeyRad = 1000*KeyCul+100*KeyNLL+10*KeyFSR+KeyISR 

100110 KeyPhy = 100000*KeyWu +10000*KeyRed+1000*KeySpn+100*KeyZ6t+10*KeyMas+KeyBra 

111 KeyTek = 100*KeySmp +10*KeyRnd+KeyWgt 

2 KeyMis = 100*KeyACC +10*Key4f +KeyMix 

— CMSENE CMS total energy 161 195 or 350 GeV 

1.16639D-5 GFERMI Fermi Constant 

128.07D0 ALFWIN alpha QED at WW treshold scale 

91.1888d0 AMAZ Z mass 

2.4974d0 GAMMZ Z width 

80.230d0 AMAW W mass 

-2.03d0 GAMMW W with. For GAMMW < RE-CALCULATED inside program 

lD-6 VVMIN Photon spectrum parameter 

0.99D0 VVMAX Photon spectrum parameter 

4 KEYDWM W- decay: 7=(ev) , 0=all ch. 

4 KEYDWP W+ decay: 7=(ev) , 0=all ch. 

-IDO WTMAX max weight for reject. ( any WTMAX<0 = is default setting) 

-1 JAKl Decay mode tau+ 

-1 JAK2 Decay mode tau- 

ITDKRC Bremsstrahlung switch in Tauola (electron muon decay) 

IFPHQT PHOTOS switch 

IFHADM Hadronisation W- 

IFHADP Hadronisation W+ 

Along with the above parameters, there are other parameters set in routine 
FILEXP; see short extract from the KORALW printout below 

* 2.03367033 W width [GeV] GAMMW 1.6 * 

* .22591156 sin(theta_W)**2 SINW2 112 * 

* Window H used only by Grace 2.0 * 

* Higgs boson parameters * 

* 1000.00000000 xpar(ll)= higgs mass amh HI * 

* 2.00000000 xpar(12)= higgs widt agh H2 * 

* DECAYS : * 

* branching ratios: * 

* .33333333 CS BR(4) IB4 * 

* .11111111 e BR(7) IB7 * 

* masses: * 

* .20000000 s AMAFIN(3) IM3 * 

* 1.30000000 c AMAFIN(4) IM4 * 
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Please note that we have used here constant Z and W widths as well as 
€\T? d^ji/ = 1 — 7j2^; KeyWu=KeyZet=l; KeyMix=2. In our approach we could 

include QCD corrections and CKM mixing as explained in ||^, vol. 1 p. 104. 
This translates simply to a change of numerical input to KORALW in setting 
the W total width and W decay branching ratios to ev and cs. We exclude, 
however, this correction from our considerations. 

2. * set 2 (modifications with respect to set 1) * 

Let us now continue with the list of corrections only. To switch spin correla- 
tions on, we leave all parameters as in previous run except: 

101110 KeyPhy = 100000*KeyWu +10000*KeyRed+1000*KeySpn+100*KeyZet+10*KeyMas+KeyBra 

3. * set 3 (modifications with respect to set 2) * 

Next, we switch on the Coulomb interaction,^ KeyRad=1000, running Z and 
W width KeyPhy=1010 and change the definition of sir? 6w The KeyMis=0 
activates a different mode of calculating sin^ 6'^ as advocated in 0] 

1000 KeyRad = 1000*KeyCul+100*KeyNLL+10*KeyFSR+KeyISR 

1010 KeyPhy = 100000*KeyWu +10000*KeyRed+1000*KeySpii+100*KeyZet+10*KeyMas+KeyBra 
KeyMls = 100*KeyACC +10*Key4f +KeyMix 

* .23103091 sin(theta_W)**2 SINW2 112 * 

Note that in this case we use for the propagators: 



S-MI + 9{s) Txs/Mx 
4. * set 4 (modifications with respect to set 2) * 



(7) 



Here we switch on all tree-level spin amplitudes, but we switch off all cor- 
rections introduced in set 3 above. In comparison to set 2, we introduce the 
following modifications: 



12 KeyMls = 100*KeyACC +10*Key4f +KeyMix 



^Note that QED coupling is taken at the low-energy limit alfinv= 137.036 
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5. * set 5 (modifications with respect to set 3) * 

Now we turn to the combined case with all tree-level spin amplitudes, running 
widths and Coulomb correction. Also the sin^ 6w is calculated by the program 
in the same way as in the case of set 3: 

10 KeyMls = 100*KeyACC +10*Key4f +KeyMix 

6. * set 6 (modifications with respect to set 5) * 

Finally with the last case we switch on the initial state bremsstrahlung and 
YFS form factor: 

1101 KeyRad = 1000*KeyCul+100*KeyNLL+10*KeyFSR+KeyISR 



Appendix B 

In this Appendix we show those input parameters used in the production of our 
Tables with the help of grc4f . They are printed on the first page of an output as an 
'echo-back', as shown below. This is for the case set 4. Description of these input 
parameters for each set will be given later. In grc4f , however, we do not prepare 
the process for the set 1. 

* 444 ########-* 

* 44444 ### ##-* 

* ###### ##— ## ##### 444—44 ##########—* 

*-####—######—##########—#### ### 444 44 ### * 

*-### ## ### ### ## 44444444444444 ### * 

*—########## ### ### 44444444444444 ### * 

*-#### ### ######—### 44 ### * 

*-#############—####### ######## 44 ######### * 

*-## ### * 

*-############ * 

* GRACE INSIDE * 

* Version 1.2 (26) * 

* Last date of change: 1996 Aug 26 * 

* Minami-Tateya Collaboration, KEK, JAPAN * 

* E*mail: grc4faiiiinajiii.kek.jp * 

* Copyright Minami-Tateya Collaboration * 

* J.Fujimoto, et.al. LEP-II Physics, WW-Generator , CERN 1996. * 



Process 


e+e 


> (3)c (4)s-bar (5)s 


(6)c-bar 


Energy : 


161. 


000 (GeV) 




Canonical CUT 




NO 




Scheme 




DN-SHELL 




sin(the_w)**2: 




0.225912 




Mass (Width) 
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* 


W-boson: 80, 


.230(2. 


03367) 


Z-boson 


: 91.189(2.49740) 


* 


u-quark : , 


.000 




d-quark 


: 0.000 


* 


c-quark: 0, 


.000 




s-quark 


: 0.000 


* 

3|c 


t -quark: 0, 


.000 




b-quark 


: 0.000 


* 


alpha : 1/128.070 


1 






* 


alpha_s: 0.12000 








* 


Experimental Cuts 










* 


Angle Cuts 










* 


Particle 


3 


-1.000 


> 


1.000 


* 


Particle 


4 


-1.000 


— > 


1.000 


* 


Particle 


5 


-1.000 


— > 


1.000 


* 


Particle 


6 


-1.000 


— > 


1.000 


* 


Energy CUT 










* 


Particle 


3 


0.000 


> 


161.000 


* 


Particle 


4 


0.000 


> 


161.000 


* 


Particle 


5 


0.000 


> 


161.000 


* 


Particle 


6 


0.000 


> 


161.000 


* 


Invariant mass CUT 








* 


Inv. Mass 3-4 


0.000 


— > 


161.000 


* 


Inv. Mass 5-6 


0.000 


> 


161.000 


* 


Inv. Mass 3-5 


0.000 


> 


161.000 


* 


Inv. Mass 4-6 


0.000 


— > 


161.000 


* 


Inv. Mass 3-6 


0.000 


> 


161.000 


* 

* 
* 


Inv. Mass 4-5 


0.000 


> 


161.000 


OPTIONS: 










* 


Calculation: 




TREE 






* 


Anomalous Coupling: 


NO 






* 


Coulomb Correction: 


NO 






* 


qCD Correction: 


NO 






* 


Gluon graph: 




NO (only for hadronic process) 


* 


Width : 




FIX 










1. * set 4 * 



We start with parameters as the set for entry 4 of our Tables and then we 
indicate respective differences only. The exact meaning of the parameters is 
described in the grc4f documentation, see Ref. 0. 



process 


= 


cCSs 


type 


= 


tree 


canon 


= 


no 


scheme 


= 


onshell 


energy 


= 


161. ODO 


alphai 


= 


128.07D0 


amw 


= 


80.23D0 


agw 


= 


2.03367D0 


amz 


= 


91.1888D0 


agz 


= 


2.4974D0 


massiv 


= 


yes 


coulomb 


= 


no 


anomal 


= 


no 


qcdcr 


= 


no 


gluon 


= 


no 


width 


= 


fixed 


opncut 


= 


O.ODO 


thecutS 


= 


180.0d0,0.0d0 


thecut4 


= 


180.0d0,0.0d0 



thecutS = 180 . OdO , . OdO 

thecute = 180. OdO, 0. OdO 

engcutS = amassl(3),w 

engcut4 = amassl(4),M 

engcutS = amassl(5),w 

engcutS = amassl(6),w 

ivmcut34= amassl(3)+amassl(4) ,w-(amassl(5)+amassl(6)) 

ivmcut56= amassl(5)+amassl(6) ,w-(amassl(3)+amassl(4)) 

ivmcut35= amassl(3)+amassl(5) ,w-(amassl(4)+amassl(6)) 

ivmcut46= amassl(4)+amassl(6) ,w-(amassl(3)+amassl(5)) 

ivmcut36= amassl(3)+amassl(6) ,w-(amassl(4)+amassl(5)) 

ivmcut45= amassl(4)+ajnassl(5) ,w-(amassl(3)+amassl(6)) 

helil = average 

heli2 = average 

heli3 = sum 

heli4 = sum 

heli5 = sum 

heli6 = sum 

end 



2. * set 5 (modifications with respect to set 4) * 
We introduce the following modifications to set 4. 



scheme = gmu 
coulomb = yes 
width = run 



The output is changed: 

* Scheme : G_mu 

* sln(the_w)**2: 0.231031 

* Coulomb Correction: YES 

* Width : RUN 



3. * set 6 (modifications with respect to set 5) * 
We introduce the following modifications to set 5. 

type = qedpsi 



The output becomes: 



Calculation: TREE + QEDPSdSR) 
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4. * set 2 * 

In its original form grc4f does not support the flag for the selection of CC-03. 
However, a simple modiflcation will enable it. We start with the parameters as 
the set for entry 4, then add an array j selg in the subroutine modmas as shown 
below, which corresponds to the case of sscc. If one sets the i-th element of 
jselg to "0", then one can omit this graph (see GRACE manual |0]). 

subroutine modmas 

implicit real*8(a-h,o-z) 

include 'incll.f 

include 'inclk.f 

* 

* 

* modification of mass and width. 



do 10 i = 1 , ngraph 
jselg(i) = 
10 continue 

jselg(57) = 1 
jselg(61) = 1 
jselg(77) = 1 

return 
end 



5. * set 3 * 



We start with the parameters as the set for entry 5 and add an array jselg 
in the same way as set 2. 
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Figure 1: The ^^ differential distribution of the "visible" ss jets where cc jets 
escape detection. The centre-of-mass energy is 161 GeV. Input parameters of type 1: 
CC-03 no spin correlation (thin line); and type 2: CC-03 spin correlations switched 
on (thick line). See Appendices A, B for a complete definition of all input parameters. 
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Figure 2: The ^^ differential distribution of the "visible" ss jets where cc jets 
escape detection. The centre-of-mass energy is 161 GeV. Input parameters of type 
2: CC-03 (thick line); and type 4: CC-43 (thin line). See Appendices A, B for a 
complete definition of all input parameters. 
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Figure 3: The ^^ differential distribution of the "visible" ss jets where cc jets 
escape detection. The centre-of-mass energy is 195 GeV. Input parameters of type 1: 
CC-03 no spin correlation (thin line); and type 2: CC-03 spin correlations switched 
on (thick line). See Appendices A, B for a complete definition of all input parameters. 
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Figure 4: The ^^ differential distribution of the "visible" ss jets where cc jets 
escape detection. The centre-of-mass energy is 195 GeV. Input parameters of type 
2: CC-03 (thick line); and type 4: CC-43 (thin line). See Appendices A, B for a 
complete definition of all input parameters. 
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Figure 5: The ^^ differential distribution of the "visible" ss jets where cc jets 
escape detection. The centre-of-mass energy is 350 GeV. Input parameters of type 1: 
CC-03 no spin correlation (thin line); and type 2: CC-03 spin correlations switched 
on (thick line). See Appendices A, B for a complete definition of all input parameters. 
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Figure 6: The ^^ differential distribution of tfie "visible" ss jets where cc jets 
escape detection. The centre-of-mass energy is 350 GeV. Input parameters of type 
2: CC-03 (thick line); and type 4: CC-43 (thin line). See Appendices A, B for a 
complete definition of all input parameters. 
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